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ABSTRACT 

We identify the optical counterpart of the ultraluminous X-ray source (ULX) NGC 1313 X-l and 
discuss constraints on its physical nature from multiband optical spectra. There is a single object 
on Hubble Space Telescope (HST) images within the aspect-corrected Chandra X-ray error circle; a 
fainter, possibly extended, feature lies near the edge of the error circle. The brighter object showed 
prominent variation in the F555W band, but was constant in the F814W band. The spectrum was 
consistent with a single power-law on 2003 Nov 17, but deviated from this on 2004 Jul 17, suggestive of 
more than one emission component. Based on the location, magnitudes, spectral shape, and variability 
of the bright object, it is likely the ULX counterpart. The red wing of the spectrum around F814W 
may be due to emission from the companion star, and the blue wing is likely from disk emission. The 
stellar population around X-l has an age older than 30 Myr, without very blue stars or young clusters. 
This places a constraint on the companion mass of the ULX as no more than 10 M Q . 
Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X- 
rays:Individual (NGC 1313 X-l) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are variable, non- 
nuclear X-ray sources found in nearby galaxies with 
isotropic luminosities higher than 3 x 10 39 erg s _1 , the 
Eddington limit of a 20Af Q black hole, roughly the 
most massive compact object that can be formed via 
core collapse of a single s tar with near solar metallic- 
ity ([Bclczvnsk i et al.ll2010f ). Growing evidence from ob- 
servations of emission line nebulae around some ULXs 
indicates that their luminosity is truly high and the 
emission is not stro ngly beamed dPakull Sz Mirionill2003l : 
IKaaret et all 120041 : IKaaret fc Corbell 120091 )" This re- 
quires accretion onto intermediate mass black holes to 
explain the high luminosity if we assume sub-Eddington 
radiation. However, the Eddington limit could be vio- 
lated under cer tain circumstances (|Watarai et al.l 120011 : 
iBegelmanl 120061 ) and mild beaming b y a factor o f a few 
is expected at high accretion rates ([King et al.1 120011 ) . 
These two factors simultaneously could explain the ma- 
jority of ULXs up to 10 41 e rg s _1 as stellar mass black 
holes (jPoutanen et al.ll2007h . The nature of ULXs is still 
uncertain. 

NGC 1313 X-l was among the first ULXs discov- 
ered with Einstein (jFabbiano fc TrinchierU Il987| ) . The 
source has shown significant spectral variability revealed 
by XMM-Newton observations. Using 12 XMM-Newton 
observ ations over a span of a few years, iFeng fc Kaaretl 
( 2006) reported a tight correlation between the X-ray lu- 
minosity and the power-law photon index of the source 
and a transition from the correlated phase to a soft state. 
This behavior is similar to that seen in some Galactic 
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black hole X-ray binaries as well as ULXs like NGC 5204 
X-l and Holmberg II X-l (|Feng fc Kaaretl [2001 . sug- 
gesting similar accretion physic s, but contrary to some 
other ULXs like NGC 1313 X-2. iDewangan etaTl (f2010h 
demonstrated that NGC 1313 X-l showed distinct tem- 
poral properties in different spectral states, manifesting a 
noisy high state and a quiet low state, interpreted as due 
to different accretion geometries. A dista nce of 4.1 Mpc 
for N GC 1313 is adopted in this paper (|Mendez et al.l 
[20021 . 

Optical observations of ULXs offer additional informa- 
tion regarding their nature, and could shed light onto 
the evolutionary history of the binary, disk geometry, and 
mode of the mass transfer (e.g.[Madhusudhan et all2008t 
IFeng fc Kaaretll2008t IKaaret fc Corbelll2009| ) . Moreover, 
a conclusion to the stellar vs. intermediate mass debate 
for ULXs would come from dynamical mass measure- 
ments, which require identification of optical counter- 
parts. NGC 1313 X-l is an important member of the 
ULX population, and has been extensively studied in X- 
rays. However, no optical counterpart of the source has 
been identified so far. Here, using data in the archive of 
the Hubble Space Telescope (HST) and Chandra X-ray 
Observatory, we report the identification of the optical 
counterpart to NGC 1313 X-l and discuss the physical 
implications. 

2. ASTROMETRY AND THE OPTICAL 
COUNTERPARTS 

NGC 1313 X-l has been observed frequently with HST 
using the Advanced Camera for Surveys (ACS) and with 
Chandra using the Advanced CCD Imaging Spectrome- 
ter (ACIS). A log of HST observations used in the paper 
is shown in Table [T] 

Chandra observation 2950, with an exposure of 20 ks, 
is the deepest in this field and thus chosen for astrome- 
try correction. The ULX is located on chip S3 of ACIS 
with a moderate offset to the optical axis of 2.4'. We 
performed source detection using the wavdetect tool in 
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TABLE 1 

HST/ACS Observations Used in the Paper 



Date 


Start Time 


Dateset 


Exposure (s) 


Filter 


2003 Nov 17 


05:52:14 a 


J8OL01010 


1160 


F814W 




06:04:50 a 


J8OL01030 b 


1160 


F555W 




07:30:30 


J8OL01040 


2520 


F435W 




09:17:54 


J8OLA1010 c 


2760 


F330W 


2004 Feb 22 


02:38:19 


J8OL05010 


2400 


F555W 


2004 Jul 17 


00:50:56 


J8PH05010 


680 


F435W 




01:08:16 


J8PH05020 


680 


F555W 




01:42:19 


J8PH05030 


676 


F814W 


2004 Oct 30 


06:57:00 


J8YY06011 


1381 


F814W 




07:25:29 


J8YY06021 


1062 


F606W 


2004 Dec 18 


21:11:19 


J8PH07020 b 


680 


F555W 



a Sub-exposures with different filters were performed alternately. 
b Observations used for astrometry correction; J8OL01030 is also 
used for photometry. 

c Observation with the aperture of the High Resolution Camera, while 
others are with the Wide Field Camera. 



CIAO 4.1.2 on a flux-corrected image created using an 
exposure map assuming a power-law spectrum with a 
photon index of 1.7 modified by Galactic absorption with 
Nh = 3.97 x 10 20 cm" 2 . We then searched in HST im- 
ages for possible optical counterparts to all X-ray sources 
detected on the S3 chip, and found one good candidate 
with both X-ray and optical emission, hereafter referred 
to as the reference object (see Figure Q] & [2]) . 

The relative astrometry between Chandra and HST 
can be improved using the reference object. Due to the 
relatively small field of view of HST, the reference ob- 
ject and X-l do not appear on any single ACS image. 
We thus created a mosaic image from two ACS obser- 
vations with the Wide Field Camera with the F555W 
filter, datasets J8PH01030 and J8PH07020. These two 
observations have an overlap of about 6 arcmin 2 . Stars in 
the overlap region were detected using the daof ind task, 
cross identified using xyxymatch, then used to calculate a 
relative shift and rotation between the two images using 
geomap. Figure [1] shows a mosaic image created using 
the multidrizzle task with these parameters. The mo- 
saic image was then alig ned to stars in the Tw o Micron 
All Sky Survey C2MASS: ISkrutskie et al.ll2006h using the 
Graphical Astronom y and Image Analysis Tool (GAIA; 
iDraper et all 12009). We find 11 2MASS counterparts 
that appear point-like, isolated, and bright (K < 15) on 
the mosaic image. Fitting to them results in a root-mean- 
square average deviation less than 2 pixels (as 0.1"), and 
an absolute astr ometric uncertainty o f 0.07" at 90% con- 
fidence (also see lFeng fc KaaretfeOOcfl . We note that this 
uncertainty does not influence the alignment of the X-ray 
and optical images. 

The reference object appears as a saturated star on 
the HST image. Its optical position is estimated from 
the cross point of the diffraction spikes as R.A. = 
03 h 18 m 18!858 decl. = -66°32'29"92 (J2000.0) with a 
negligible statistic error of about 0.01". It has 82 pho- 
tons detected on the Chandra image with a significance 
of 38cr and an error radius of 0.17" at 90% confidence. 
Using MARX simulation, the systematic offset in the po- 
sition calculated by wavdetect caused by asymmetry of 
the photon distribution is estimated as AR.A. = 0.03" 
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Fig. 1. — The mosaic HST image in F555W containing X-l and 
the reference object used for astrometry correction. The arrow 
points north and has a length of 1'. 




Fig. 2. — Optical image of the reference object. The plus and the 
circle indicate the Chandra position and error circle of 0.6" radius, 
respectively. The cross indicates the optical center of the source. 
The arrow points north and has a length of 1". 



and Adecl. = —0.04". These values were used to correct 
the offset in the coordinates from wavdetect. In Fig- 
ure [21 the optical image of the reference object is shown. 
The cross indicates the optical center and the plus indi- 
cates the X-ray position with the circle showing the 0.6" 
radius absolute astrometry uncertainty of Chandra. The 
X-ray image can be registered onto the optical image by 
aligning the plus and the cross. 

For X-l, there are 2968 photons detected with Chan- 
dra, and the statistic error radius is 0.04" at 90%. The 
systematic errors are AR.A. = —0.08" and Adecl. = 
0.02" and are taken into account. The corrected po- 
sition of X-l relative to HST is R.A. = 03 h 18 m 20?013 
decl. = -66°29T0"79 (J2000.0). The uncertainty on this 
corrected position is similar to the statistical error on the 
X-ray location of the reference object, i.e. 0.17", since all 
other errors are significantly smaller. However, the ro- 
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Fig. 3. — HST images around NGC 1313 X-l in different bands. On the F555W image the dashed and solid circles indicate, respectively, 
the original and corrected Chandra position errors of X-l. On F435W, the circle indicates a group of stars nearby. The horizontal bars 
point to the bright optical object in the corrected X-ray error circle. The F606W image is from the observation on 2004 Oct 30, and others 
are obtained from observations on 2003 Nov 17. All arrows point north and have a length of 1". 



tation of the X-ray image cannot be fixed with a single 
object, and the accuracy in rotation of Chandra astrom- 
etry is not discussed in any available Chandra documen- 
tation. We estimated the rotation accuracy from a few 
Chandra observations with multiple optically identified 
X-ray sources and obtained a typical rotation about 2' 
between Chandra and 2MASS. Taking this value into ac- 
count, with a distance of ^200" between the reference ob- 
ject and X-l, the error caused by rotation is about 0.1". 
Combining the two errors, 0.1" and 0.17", in quadra- 
ture, we estimate that the corrected X-ray position has 
an uncertainty of about 0.2" at the 90% confidence level. 

The F555W image around X-l is shown in Figure 02 on 
which a large, dashed circle indicates the original Chan- 
dra position with an absolute error of 0.6", while a small, 
solid circle indicates the corrected relative position with 
a radius of 0.2". The astrometry correction successfully 
rules out a few objects as possible counterparts to X-l, 
but leaves one relatively bright object within the cor- 
rected relative error circle and another faint feature just 
to the east of the circle. Even if the error radius of the 
corrected Chandra position is 2.5 times larger, i.e. 0.5", 
this conclusion does not change. All available HST im- 
ages of X-l in other bands (F330W, F435W, F606W, 
and F814W) are also displayed with the same scale as 
the F555W image. Horizontal bars are plotted pointing 
to the bright object in the error circle. The bright ob- 
ject is obvious in all bands. The faint feature is more 



significant at shorter wavelength. 

At about 2" to the southwest of X-l, a group of stars 
with an extent of about 1.5" (30 pc) are obviously seen, 
indicated by a circle on the F435W image. A few bright 
objects in the group may consist of multiple unresolved 
sources, especially on images of long wavelength. 

3. PHOTOMETRY 

3.1. Flux and Variability of the counterpart 

The region around X-l on the HST images is crowded, 
so we decided to perform aperture photometry with an 
aperture radius of 0.1", in order to reduce contamina- 
tion from nearby sources. To test the aperture correc- 
tion for this relatively small aperture, we selected nearby 
isolated bright objects and performed photometry us- 
ing both small and large apertures. We found a typical 
uncertainty around 0.02 mag, which is added into the 
uncertainty in our reported magnitudes. We note that 
PSF photometry in this circumstance does not offer more 
reliable measurement, as features near the source have 
somewhat complex morphology and are not point-like. 
Background levels were estimated from a concentric an- 
nulus where no sources appear; if a couple dim sources 
could not be avoided, we then chose a background re- 
gion as large as possible to depress the source contri- 
bution. The intrinsic source flux or magnitude was cal- 
culated using the SYNPHOT package using the count 
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TABLE 2 

Magnitudes corrected for Galactic extinction of the 
bright optical object in the corrected error circle of x-l. 



Band 2003 Nov 17 2004 Feb 22 2004 Jul 17 2004 Oct 30 



F330W 22.26 ±0.03 

F435W 23.00 ±0.03 ■■• 23.11 ±0.04 ••■ 

F555W 23.65 ±0.04 23.68 ± 0.04 23.90 ± 0.05 

F606W ■ ■ ■ ■ ■ • ■ • • 24.08 ± 0.04 

F814W 24.92 ±0.06 ■■• 24.86 ± 0.06 24.92 ±0.07 



Vega magnitudes in Johnson-Cousins bands 



U 22.62 ±0.03 

B 23.64 ±0.03 ■■• 23.74 ±0.04 

V 23.72 ±0.04 23.74 ±0.04 23.96 ± 0.05 23.85 ± 0.04 

I 23.65 ±0.06 ■■■ 23.59 ±0.06 23.65 ±0.07 




14.5 14.6 14.7 14.8 14.9 14.5 14.6 14.7 14.8 14.9 
logf 

Fig. 4. — Multiband spectrum of the bright optical object in 
the corrected error circle of X-l from different epochs. The lines 
in different panels are the same, indicating the best-fit power-law 
model, F v <x y0.8ii0.07 to the data on 2003 Nov 17. 
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Fig. 5. — Light curves in the F555W band for the bright opti- 
cal object in the corrected error circle of X-l and another 8 nearby 
sources within 10" with similar magnitude. The instrumental mag- 
nitude is computed as m = 25 — 2.5 log r, where r is the observed 
net count rate. The top one is the counterpart of X-l, and the 
magnitudes for other sources are shifted vertically for clarity. The 
dashed lines are the best fit constants for each. 

Figure O All these 8 sources show a light curve con- 
sistent with a constant, with \ 2 = 0.3 — 3.3 for 2 d.o.f. 
Also, the multiband spectra of these 8 sources are closely 
consistent between the two observations on 2003 Nov 17 
and 2004 Jul 17. Therefore, we conclude that the varia- 
tion in the F555W band for the counterpart is intrinsic. 
We also checked the variability in F435W and F814W 
for the ULX and comparison sources in the same way. 
Within the photometric errors, there is no evidence for 
variability of the ULX. However, the data in the F435W 
band are incompatible with the level of variability seen 
in the F555W band only at the 90% confidence level. 
For the faint object to the east of the corrected er- 



rate obtai ned from photometr y corrected to an infinite 
aperture (|Sirianni et al.l I2005D and Ga lactic reddening 
E[B - V) =0.11 (|Schlegel et al.l (19981) . The PSF may 
vary slightly with temperatur^l, causing a systematic 
shift about 0.02-0.05 mag between observations with such 
a small aperture. Several nearby isolated bright stars 
were used to calculate a zero point shift to correct for 
this effect using the magnitudes measured on 2003 Nov 
17 as reference. The extinction corrected fluxes at differ- 
ent bands were fitted to a power-law spectrum, F v oc v a , 
where v is the pivot frequency of the filter in Hz and F v 
is the flux in erg s _1 cm~ 2 Hz -1 . The same power-law 
spectrum was assumed during the conversion from rate 
to flux; a couple cycles of iteration were needed here to 
produce a consistent power-law index before and after 
the conversion. 

ST magnitudes corrected for Galactic extinction for 
the bright object are listed in Table [2] They are also 
transformed into Vega magn itudes in the closest Johnson 
UBV and Cousins I bands ([Sirianni et al.ll2005f ). Addi- 
tional errors introduced by the uncertainty on the as- 
sumed power-law index is taken into account but are not 
significant. The source has an absolute visual magnitude 
M v = -4.35 ± 0.04 and color (B - V) = -0.08 ± 0.05 
(Johnson magnitudes with VEGA zero points corrected 
for Galactic extinction) from observations on 2003 Nov 
17. The intrinsic fluxes in different bands from differ- 
ent epochs are plotted in Figure 0] where the solid lines 
indicate the best power-law model fitted to the spec- 
trum on 2003 Nov 17. As shown in Figure [4] the fluxes 
measured on 2004 Jul 17 and Oct 30 show a devia- 
tion from the best-fit power-law obtained from 2003 Nov 
17, except in the F814W band. The largest variation 
is seen in the F555W filter with a magnitude increase, 
Am = 0.25 ± 0.06 (here the aperture correction error 
is not included as we are considering relative shifts for 
the same filter), for the observation on 2004 Jul 17 with 
respect to that on 2003 Nov 17. Fitting to a constant 
with the three F555W magnitudes results in \ 2 = 23.4 
for 2 degrees of freedom (d.o.f.). To investigate whether 
the variation is intrinsic or instrumental, we selected 8 
sources with similar magnitude in F555W within 10" of 
X-l, and checked their magnitude variation (again, with- 
out considering the aperture correction error), shown in 

5 http:/ /www. stsci.edu/hst/acs/documents/isrs/isr0712.pdf 
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ror circle, it has a reddening corrected ST magnitude 
of 25.3 ± 0.1 in F555W on 2003 Nov 17. Interestingly, 
there seems to be an extended feature connecting it to 
the bright object, as one can see at long wavelengths, in 
particular in F814W, while the object itself, at the tip of 
this extended feature, becomes invisible in F814W. This 
faint object shows a constant flux within errors. Due to 
the deepness of HST observations, such a faint object 
has a relatively high chance probability to show up in 
any circle of 0.2"-radius. We currently cannot constrain 
its nature or rule out the possibility that it is not the 
physical counterpart to X-l. However, due to two facts 
we prefer the bright object as the better candidate of 
the counterpart to X-l. First, its magnitude and spec- 
tral shape are similar to the optical counterparts of most 
other ULXs tha t have a unique identification in optical 
(|Tao et al.ll2011| ). Also, variability is expected for an ac- 
creting source. 

3.2. X-ray to Optical Flux Ratio 

Two sets of HST observations used here have simulta- 
neous Chandra observations. Unfortunately, X-l caused 
severe pileup on the CCD with about 25% of events be- 
ing piled from multiple photons in these observations, 
hampering us from obtaining reliable spectra and fluxes. 
Based on the recorded X-ray count rates and using 
PIMMS, we find that the X-ray flux of the source dur- 
ing the two Chandra ob servations is within t he range 
found by XMM-Newton (|Feng fe Kaaretl [2001 . There- 
fore, we adopt this range as the X-ray flux to calculate 
the X-ray to optical flux ratio, defined as log(/x//v) = 
log/x + mv/2.5 + 5.37, where fx is the observed X- 
ray flux in the 0.3-3.5 keV band and my is the ob- 
served visual magnitude (jMaccacaro et al.lll988| ). Insert- 
ing fx = (1.4 -4.6) xl0~ 12 er g 
flux range in lFeng fe Kaaretl 

have log(/x/ fv) — 3.1 — 3.7. This is significantly larger 
than the X-ray to optical flux ratio of active galactic nu- 
clei ( A GN) , which is usually found in the range f rom —1 
to 1.7 (jMaccacaro et al.lll988tlStocke et al.lll99lD . 

For the purpose of distinguishing low mass and high 
mass X-ray binaries, an alternative definition of the ra- 
tio, called X-ray to optical color index, has been proposed 
as £ = Bo + 2. 5 log Fx, where B is the extinction cor- 
rected B magnitude and Fx is the observed X-ray flux 
density in 2-10 keV in units of /j Jy (jvan Paradiis 1981; 
Ivan Paradiis fe McClintockl IT995h ■ With ffp = 23.6 and 
Fx = 0.03 — 0.23 (LtJy (also derived from lFeng fe Kaaretl 
2006), £ is estimated to be in the range of 20-22 for NGC 
1313 X-l. For high mass X-ray binaries during outbursts, 
£ is around 12-18; for low mass X-ray binaries (LMXBs), 
the £ d istribution peaks around 21-22, w ith a deviation 
of 1-2 (|van Paradijs fe McClintockl 119951 ). Therefore, in 
terms of £, the ULX is more similar to LMXBs. 

3.3. Color- Magnitude Diagram 

The magnitudes for objects within 10" of X-l (about 
200 pc correspondingly) were calculated using the 
DAOPHOT package in the Interactive Data Language 
(IDL). A total number of 123 bright sources with a flux 
4cr above the background were detected using daof ind 
from observations on 2003 Nov 17 in the F435W and 
F555W bands. All these sources have a local fit x 2 < 3 



r g cm " s x , found from the 
I (20061 ). and m v = 24.1, we 
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Fig. 6. — Color-magnitude diagram for objects within 10" 
(200 pc) of NGC 1313 X-l. Dashed lines are theoretical isochrones 
at ages of 10 7 ' 2 , 10 7 ' 5 , and 10 7 ' 8 yr, respectively. Filled circles are 
for objects in the nearby group of stars, i.e. inside the circle on 
F435W in Figured] 



to exclude extended sources and image defects. We 
performed PSF photometry for these sources and cre- 
ated a color-magnitude diagram shown in Figure [H for 
F555W vs. F435W— F555W. The magnitudes are com- 
puted assuming an AO spectrum based on their mean 
color, normalized to VEGA zero point, and corrected 
for the Galactic extinction. Typical errors are 0.03 for 
magnitudes and 0.05 for colors. For the nearby group 
of stars, eight bright, point-like sources are detected and 
displayed in the figure as filled circles. 

Observations of H n regions in NGC 1313 suggest a flat 
distribution of oxyg en abundance with 12 + log(0/H) = 
8.4 across the disk (IWalsh fe Rovlfl997l) . converting to a 
metall icity of Z = 0.006 based on standard solar compo- 
sition (|Grevesse fe S auval 19981). Theoretical isochrone^ 
of stars with the same metallicity at ages o f 10 72 , 10 75 , 
and 10 7 - 8 yr are plotted (jMarigo et al.ll2008[ ). As one can 
see, most objects are older than 10 7 5 yr, and almost none 
of them is younger than 10 7 2 yr. In the diagram, stars 
in the nearby group do not appear differently from field 
sources. We also checked the diagram between F555W 
and F814W and obtained the same conclusion. 

4. DISCUSSION 

We have demonstrated the identification of the opti- 
cal counterpart of NGC 1313 X-l as the brighter object 
in the corrected error circle. Deviation from a single 
power-law spectrum on 2004 Jul 17 may suggest that 
at least two components are responsible for the optical 
emission. Due to a constant flux in F814W and variabil- 
ity in F555W, the red wing of the spectrum is possibly 
dominated by emission from the companion star and the 
blue wing is likely dominated by disk radiation. 

Assuming that the red component is mainly due to star 
light from the companion, we can place constraint on its 
absolute magnitudes My > —4.1, Mj = —4.5, and color 
(V — I)o > 0.4, where the values are translated from 
F555W and F814W on 2004 Jul 17 and the inequality 
accounts for the fact that the F555W flux of the red 
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component is lower than the observed value. A K5-M0 
II bright giant would match the /-band magnitude and 
be consistent with the color. The F814W flux on 2003 
Nov 17 should also be due to the companion emission 
under this hypothesis, and the single power-law is then 
a composition of two components. The distinct optical 
spectra between 2003 Nov 17 and 2004 Jul 17 indicate 
that the disk emission has a varying break frequency w 
(5 — 8) x 10 14 Hz, corresponding to a break temperature 
lb = hvh/k ~ 3 x 10 4 K, where h and k are the Planck 
and Boltzmann constants, respectively. 

Disk emission in the optical can be decomposed into 
two origins, one from intrinsic viscous release of the 
multicolor disk, and the other from X-ray reprocessing 
(irradiation) in the outer regions of the disk. The in- 
trinsic emission of a standard thin disk has a spectral 
form of f 1 / 3 at frequencies above kT out /h, where T out is 
the disk temperature at the outermost radius, and then 
breaks to v 2 at much lower frequencies. The two point 
power-law index between the F555W and F435W fluxes 
on 2004 Jul 17 is a = 1.4 ± 0.3 corrected for Galactic 
extinction; if additional extinction exists in NGC 1313, 
a would be even larger. Thus, the blue wing of the 
spectrum is steeper than v 1 ^ and cannot be explained 
by intrinsic disk emission only, unless the disk is trun- 
cated with T out ~ It,. For viscous dissipation on a stan- 
dard thin disk, we have R ou t/Rin = (Tout /Tin) " 4 / 3 - If 
we assume T in = 0.1-1 kcV (10 6 - 10 7 K) for either 
intermediate or stellar mass black holes, then we get 
Ront/Rm = 100 — 2000, suggesting the presence of a rel- 
ative ly narrow disk, as h as been propo sed for NGC 6946 
X-l JKaaret et alJl20tnh and SS 433 (|Perez fc Blundelll 
I2010T ). Otherwise, if disk irradiation is dominant, it is 
expected to have an effective temperature Tj r » Tb with 
a surface area around 10 24 cm 2 to reach the observed flux 
density. Therefore, a disk size about 10 12 cm is inferred. 

As optical fluxes at different bands were not obtained 
simultaneously, an alternative explanation of the non- 
power-law spectrum seen on 2004 Jul 17 could be due to 
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short-term variation of the disk emission. The last sub- 
exposure of F435W started in advance of the F555W 
observation by 544 s, which, corresponding to a light 
cross distance over 10 13 cm, is sufficiently long to re- 
spond to the rapid X-ray variation. Strong, fast optical 
variations have been found for both LMXBs and ULXs 
(|Bradt fc McClintockjH98l IGrise et al1l200ll . and NGC 
1313 X-l is highly vari able in X-rays at short timescales 
(|Dewangan et al.ll2010ft . Simultaneous X-ray and optical 
observations in the future could test different scenarios. 

The properties of the companion star could be con- 
strained by comparing the observed My and (B — V)q 
with theoretical values derived from binary evolution 
(jPatruno fc Zampieri|[2008l[20Tl . The data can be best- 
fitted with a donor of 10-15 Mq during hydrogen shell 
burning given a black hole mass of 10-100 Mq, or dur- 
ing main sequence only if the black hole mass is close to 
100 Mq. Also, it can be ruled out that the donor starts 
to fill the Roche-lobe during the gi ant phase if the black 
hole mass is greater than 20 Mq (jPatruno fc Zampieril 
l2010h . 

The color-magnitude diagram and isochrones reveal 
a stellar population around NGC 1313 X-l with most 
stellar ages greater than 10 75 Myr. The nearby region 
(within 10" or 200 pc) also lacks very blue stars and 
young star clusters. This places a constraint on the mass 
of the companion star as no more than ~ 10M Q . This is 
consistent with the hypothesis above that the companion 
may be a late type bright giant. The ULX is close to a 
group of stars, which has an extent of about 30 pc and 
has a similar age to field stars. Whether or not the ULX 
belongs to the star group is unknown. This is different 
from the environment around NGC 1313 X-2, which is 
found to be associate d with a young gr oup consisting of 
bright and blue stars (jGrise et al.ll2008h . 

We thank the referee for useful comments that have 
improved the paper. 

Facility: CXO(ACIS), HST(ACS) 



Marigo, P., Girardi, L., Bressan, A., Groenewegen, M. A. T., 
Silva, L., & Granato, G. L. 2008, A&A, 482, 883 

Mendez, B., Davis, M., Moustakas, J., Newman, J., Madorc, B. 
F., Freedman, W. L. 2002, AJ, 124, 213 

Pakull, M. W., & Mirioni, L. 2003, Rev. Mex. AA Ser. Conf., 15, 
197 

Patruno, A., & Zampieri, L. 2008, MNRAS, 386, 543 

Patruno, A., & Zampieri, L. 2010, MNRAS, 403, L69 
Perez, S. M., & Blundell, K. M. 2010, MNRAS, 408, 2 
Poutanen, J., Lipunova, G., Fabrika, S., Butkevich, A. G., & 

Abolmasov, P. 2007, MNRAS, 377, 1187 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Schlegel, D. J., Finkbeiner, D. P., Davis, M. 1998, ApJ, 500, 525 
Sirianni, M., et al. 2005, PASP, 117, 1049 

Stocke, J. T., Morris, S. L., Gioia, I. M., Maccacaro, T., Schild, 
R.. Wolter, A., Fleming, T. A., & Henry, J. P. 1991, ApJS, 76, 
813 

Tao, L., Feng, H., Kaaret, P, & Grise, F. 2011, ApJ, submitted 
van Paradijs, J., & McClintock, J. E. 1995, in X-Ray Binaries, ed. 

W. H. G. Lewin, J. van Paradijs, & E. P. J. van den Heuvel 

(Cambridge: Cambridge Univ. Press), 58 
van Paradijs, J. 1981, A&A, 103, 140 
Walsh, J. R., & Roy, J. R. 1997, mnras, 288, 726 
Watarai, K.-Y., Mizuno, T., & Mineshige, S. 2001, ApJ, 549, L77 



